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Abstract 

The problem of computing the thermodynamic properties of a one- dimensional gas 
of particles which transform in the adjoint representation of the gauge group and 
interact through non-Abelian electric fields is formulated and solved in the large N 
limit. The explicit solution exhibits a first order confinement-deconfinement phase 
transition with computable properties and describes two dimensional adjoint QCD 
in the limit where matter field masses are large. 
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Two dimensional quantum chromodynamics (QCD) with adjoint representation 
matter fields is the simplest field theoretical model which exhibits some of the 
important common features of string theory and the confining phases of gauge 
theory. Most notable are the infinite number of asymptotically linear Regge 
trajectories and a density of states which increases exponentially with energy. In 
two dimensions, the Yang-Mills field itself has no propagating degrees of freedom. In 
adjoint QCD, the matter fields provide dynamics by playing a role analogous to the 
transverse gluons of higher dimensional gauge theory. In fact, dimensional reduction 
of three dimensional Yang Mills theory produces two dimensional QCD with massless 
adjoint scalar quarks. Moreover, since adjoint matter fields do not decouple in the 
infinite iV limit, the large N expansion is of a similar level of complexity to that of 
higher dimensional Yang-Mills theory. One would expect it to exhibit some of the 
stringy features of the confining phase which are emphasized in that limit. 

Although adjoint QCD is not explicitly solvable, even at infinite N, details 
of its spectrum were readily analyzed by approximate and numercial techniques 
HI 0, 0, ||]. In addition, Kutasov exploited an argument which was originially 
due to Polchinski || to show that the confining phase must be unstable at 
high temperature and suggested it as a tractable model where the confinement- 
deconfinement transition could be investigated. 

In this letter, we shall formulate and find an explicit solution of the large N 
limit of a simplfied version of adjoint QCD. We shall consider a one-dimensional 
gas of non-dynamical particles which have adjoint color charges and which interact 
with each other through non-Abelian electric fields. Because there are no dynamical 
gluons which could screen adjoint charges in one dimension, at low temperature and 
density, adjoint quarks are confined^. They form colorless "hadron" bound states 
with two or more adjoint quarks connected by non-dynamical strings of electric flux. 
The large N limit resembles a non-interacting string theory in that, at infinite N, the 
energy of a state is proportional to the total length of all strings of electric flux plus 
a chemical potential times the total number of quarks. The property of confinement 
is defined by estimating the energy required to introduce an external fundamental 
representation quark-antiquark pair into the system. In the confining phase, where 
the hadron gas is dilute, the quark-antiquark energy is proportional to the length 
of the electric flux string which, to obtain gauge invariance, must connect them. 
This gives them a confining interaction. Some typical confined configurations are 
depicted in Fig. 1. In the confined phase, the average particle number density and 
the energy density are small — in the large N limit both are of order one, rather 
than iV 2 which one would expect from naive counting of the degrees of freedom. 
This is consistent with the fact that in a confining phase the number of degrees of 
freedom, i.e. hadrons, is independent of N. In contrast, in the deconfined phase, 
since the number of degrees of freedom, i.e. quarks and gluons, is proportional to 
TV 2 the particle density and energy are also of order N 2 . 

As temperature or density is increased, eventually we arrive at the situation 
where there are electric flux strings almost everywhere. Then, adding an additional 

^dn higher dimensions, an adjoint charge and a gluon could form a color singlet bound state. 
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Figure 1: Some examples of states in the confining phase. Each adjoint quark 
connects with two electric flux strings and fundamental quarks in (b) connect with 
one string. 
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Figure 2: Typical configurations in the deconfined phase. 



flux string, or modifying the existing network of strings to accommodate a funda- 
mental representation quark-antiquark pair, involves a negligibly small addition to 
the energy of the total system (see Fig. 2). This is typical of the deconfined phase. 

Between these two phases is a transition, which we shall show in this paper, 
is of first order. In the string picture, this phase transition occurs when the 
strings in a typical configuration percolate in the one- dimensional space. The order 
parameter is the Polyakov loop operator || [7| which measures the exponential of the 
negative of the free energy which is required to insert a single, unpaired fundamental 
representation quark source into the system. This free energy is infinite (and the 
expectation value of the Polyakov loop is zero) in the confining phase and it is finite 
in the deconfined phase. 

In the Hamiltonian formulation of two dimensional Yang-Mills theory, the 
electric field is the canonical conjugate of the spatial component of the gauge field, 
A a (x),E b (y) 



iS ab (x — y). The Hamiltonian isQ 

2 N 2 

H = jdx e -Y,(E\x)f 



a=l 



and the Gauss' law constraint, which takes the form of a physical state condition, is 

K 



( d K 

-j-E a (x) - f ahc A\x)E c (x) + Y. T i 8( y 
\dx i=l 



X — Xi)\ ^ 
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(2) 



2 Here, for concreteness, we consider U(N) gauge theory. The gauge field A — A a t a , with t a the 
generators in the fundamental representation. 
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There are particles with adjoint representation color charges located at positions 
Xi, . . . ,xk- T? are generators in the adjoint representation operating on the color 
degrees of freedom of the fth particle. 

In the functional Schrodinger picture, the states are functionals of the gauge field 
and the electric field is the functional derivative operator E a (x) = \ SA t^ ■ The 
functional Schrodinger equation is that of a free particle 

r ( e 2 N ' 2 5 2 \ 

J dX [~2 S W [A; Xl > ■ ■ ■ ' Xk] = 8 [A;Xl > ■ ■ ■ ' Xk] 

(3) 

Gauss' law implies that the physical states, i.e. those which obey the gauge 
constraint (|2]), transform as 

[A"; x u ...,x K }= g&fa) . . . 9% bK M* bl ~ bK [*> *i> • • • . x *\ ( 4 ) 

where A 9 = gAg^ + igVg^ is the gauge transform of A. For a fixed number of 
external charges, this model is explicitly solvable. In the following we shall examine 
its thermodynamic featurs, where we assume that the particles have Maxwell- 
Boltzmann statistics. 

We find it convenient to work with the grand canonical ensemble. The partition 
function is constructed by taking the trace of the Gibbs density e~ H l T over physical 
states. This can be implemented by considering eigenstates of A a (x) (and an 
appropriate basis for the non-dynamical particles) |-A)e ffll . . . e aic . Projection onto 
gauge invariant states involves a projection operator which has the net effect of 
gauge transforming the state field at one side of the trace, and integrating over all 
gauge transformations ||. The resulting partition function is 

Z[ Xi ,T\ = J[dA][dg] (A\e- H / T \A°)Trg Ad (x 1 )...Tr g Ad (x K ) , (5) 

where [dg(x)] is the Haar measure on the space of mappings from the line to 
the group manifold and [dA] is a measure on the convex Euclidean space of 
gauge field configurations. For U(N), the trace in the adjoint representation is 
Tr g Ad (x) = | Tr g(x) \ 2 where g{x) is in the fundamental representation. In order 
to form the grand canonical ensemble, we average over the particle positions by 
integrating / dx\ . . . J dxx, multiply by the fugacity to the power K, X K , divide by 
the statistics factor 1/K\ and sum over K. The result is 

Z[X,T] = J[dA}[dg} e ~ s ^ [A ' 9] (6) 

where the effective action is 

e S cS [A,g] = ^| e -H/T \ A gj exp (J & x \ |Tr gf^j (7) 

The Hamiltonian is the Laplacian on the space of gauge fields. Using the explicite 
form of the heat kernel 

(A\ e~ H/T \A 9 } ~ exp (- J dx J Tr (A - A 9 ) 2 ^ , 
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we see that the effective theory is the gauged principal chiral model with a quadratic 
potential 

S#[A,g) = Jdx (^Tr|V£ + z[As]| 2 -A|Tr £| 2 ) (8) 

This effective action with A = was discussed by Grignani et.al. [[|. It 
is gauge invariant, S e s[A, g] = S e g[A h , hgh^}, and has the global symmetry 
SeffLA, g] = S e g[A, z g], where z is a constant element from the center of the 
gauge group, which for U(N) is U(l) and for SU(N) is Zjy. 

The realization of this center symmetry governs confinement || |7|]. When 
the symmetry is represented faithfully, the theory is in the confining phase. The 
Polyakov loop operator Tr g(x) transforms under the center as Tr g(x) — > z Tr g(x). 
Thus the expectation value of the Polyakov loop operator must average to zero if 
the symmetry is not spontaneously broken. This expectation value is interpreted as 
the free energy of the system with an additional external charge in the fundamental 
representation of the gauge group situated at point x, F[x, A,T] = —Tin (Tr g(x)). 
For finite N, and D = 3, ideas of universality have been applied to study phase 
transitions with this order parameter in SU(N) gauge theory [ ITS | . The phase 



transition should be second order for N = 2 and first order for N > 2. 

Here we amalyze the effective theory (|6]),(|8|) in the large N limit. If we rescale the 
coupling constant so that %r = jf, both terms in the action @ are of order N 2 and 
in the large N limit the partition function is dominated by the configuration which 
minimizes the action. Gauge invariance can be used to diagonalize the matrices 
gij(x) = e loli ^5ij. The density of eigenvalues p(9,x) = ^ X^i — cti(x)) 
corresponding to the large N saddlepoint now characterizes the properties of the 
system. A constant density p coa f,(9, %) = ^~ realizes the center symmetry, and thus 
corresponds to the confining phase. A density peaked at some value of 9 explicitely 
breaks the center symmetry, and corresponds to a deconfined phase. 

If, in the general case, we consider the Fourier expansion 

p(9,x) = ^ + ^J2^(x)e- md , c n (xY = c_ B (x) , (9) 
Ztt Ztt ,„ 

the Fourier coefficients c n (x) characterize the possible deconfined phases of the 
theory. If one of them were non-zero, we would have in the infinite iV limit 
j. ^ r g n (x)) = c n (x). This would indicate that a composite of n fundamental 
quarks would have finite free energy and would not be confined. 

In order to find the configurations of the eigenvalue density (§) that minimize 
the action, we shall use the collective field theory approach of Refs. [[II], 0, |T3|1 . 



Alternatively to the gauge fixing that we have discussed, we consider (§) in the gauge 
A = which can be fixed on the open line. Then the thermodynamic problem is 
equivalent to unitary matrix quantum mechanics 



Z[X, T] = J [dg] exp (- J dx f^Tr | Vg\ 2 - A |Tr g\ 



(10) 



This model can be solved in the large iV limit by the methods of collective field 
theory. The method is essentially based on the relation between matrix quantum 
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mechanics and nonrelativistic fermions p^| . In the large N limit the eigenvalue 
density obeys a classical, saddle point equation which can be deduced from canonical 



analysis of the collective field theory Hamiltonian 1TTJ 



H 



<:W 
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and subsequent Wick rotation to imaginary time. Here n(6 l ) is the variable which is 
the canonical conjugate of p(9), so that the Poisson bracket is {p{9), 11(0')} = 5(9 — 
9'). The velocity of the Fermi fluid is v(9) = 811/89. In the equations of motion 
following from (|TT[) , we change t — > ix, v — > —iv and obtain 



dp , 8 

Tx +1 m {pv) 







(12) 



8v 



dv 



8x + 1V ~89 



7T- 7 p^ + 2AIm(e- l %(x)j =0. (13) 

It is expected that the solution of these equations corresponding to an equilibrium 
state of the system is a constant p(x,9) = po(9). At least at sufficiently low 
temperature or, equivalently, at sufficiently large 7, the system is in the confining 
phase with unbroken center symmetry, so that p = p con f. = 1/27T. This is always a 
solution of the equations of motion (|T2| ) ,([T3| ) since c\ = 0. 

However, this solution is stable against small fluctuations only if 7 is large 
enough. To find the spectrum of excitations in this phase, we linearize the equations 
of motion around p con f.- To do this, we consider the c n (x) of eq. (^) and v 
infinitesimal. The resulting equation for c n (x) is 



7 2 rz 2 



n,— 1 , 



c n (x) =0 ; n^O 



(14) 



At 7 = 7c(A) = 4A, the lowest eigenvalue corresponding to n — ±1 goes to zero. For 
smaller 7 this eigenvalue is negative and the strong coupling solution is unstable with 
c±i the first modes to become unstable. However, for reasons which will become 
clear once we consider the weak coupling phase, 7c (A) should not be identified with 
the point of the deconfining phase transition. 

The solution in the deconfined phase can be obtained by integration of eq. ( |T3| ) 
at v = 0. The density po(0) can always be chosen to be an even function of 9. Thus 
C\ is real, and one finds from eq. (|T2"D: 



I 2 

pJQ) = - -E + 2\ Cl cos9 . 
ix y 7 



Outside the region [— 6* max , 9 max ] with # max = ir — arccos -j^, the density p(9) is zero. 

The Fermi energy E and the constant c\ are to be determined from the 
normalization condition /_ e max dO po(9; E,ci) = 1, and the consistency condition 
/-<? dO cos 9 po(9; E,ci) = c\ derived from eq. (Q). It follows from these equations 
that # max tends to zero at 7 — >• and grows with the increase of 7. Eventually 
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Figure 3: The large N phase 
coupling (confining) phase, II 
which the weak coupling phase 
phase transition: 7 = 70(A), 3 - 
7 = 7c(A) 



diagram of the one- dimensional model. I - strong 
- weak coupling (deconfining) phase; 1 — line on 
terminates: 7 = 7* (A), 2 — line of the first-order 
- line of the instability of the strong coupling phase: 



it reaches ir, where the weak coupling phase terminates, because the eigenvalue 
distribution begins to overlap with itself due to 2-7r-periodicity. At the critical point 
E* = 2Aou, the normalization and consistency integrals can be done explicitly. We 
find that c u = 1/3 and 7*(A) = |f| A « 4.324 A. 

We obtain the following picture of the deconfining phase transition (Fig.3). The 
weak and strong coupling phases can coexist, because 7 C (A) < 7* (A), although 
the region, where both phases are stable is very narrow, since 7 C (A) and 7* (A) are 
numerically close to each other. The phase transition is of the first order and takes 
place at some 70(A) between 7 C (A) and 7* (A). At the point of the phase transition 
the free energies of both phases are equal to each other. Substituting po(6) into 
equation ( jTT|) one can find the free energy per unit volume, to leading order in the 
large N limit, 

F ( 0, in the confining phase , , 

N 2 I \E — \\c\ — 24, in the deconfined phase 

The equations determining the critical line can be solved numerically to obtain 
70(A) = 4.219 A. 

The model which we have considered in this Section is adjoint QCD in the limit 
where the particles are heavy. The fugacity parameter can be computed from a 
one-loop diagram as A = yj ^e _m//T , the exponential being simply the Boltzmann 
weight of a particle with mass m. It is assumed that m» T for classical statistical 
mechancis to be applicable and m >> e to suppress pair production. Our results 
indicate that the phase transition is of first order with critical line approximately 
given by the equation « 4.2-y/^e _m//T . There exists a region of parameters 
in which this equation has a solution and the conditions of applicability of our 
simplified model are satisfied. 

The work of K. Zarembo was supported in part by INTAS grant 94-0840. The work 
of G. Semenoff and O. Tirkkonen was supported in part by NSERC of Canada. 
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